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Abstract 

We present excited states density functional theory (DFT) to calculate band gap 
for semiconductors and insulators. For the excited states exchange-correlation func- 
tional, we use a simple local density approximation (LDA) like functional and it 
gives the result which is very closed to experimental results. The linear muffin-tin 
potential is used to solve the self consistent Kohn-Sham equation. 
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1 Introduction 

First-principles calculations based on the density functional theory (DFT) 
have been eminently successful for the study of ground state properties of 
electrons in a solid [1]. However, it has always been understood that the 
spectrum and wave-functions of the associated one-electron like Kohn-Sham 
equation have no specific significance beyond the fact that they are used to 
obtain the ground-state density [2]. Thus the unoccupied part of the Kohn- 
Sham spectrum cannot be expected to describe correctly the excitations of a 
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many-electron system. On the other hand, the energy eigenvalue of the high- 
est occupied orbital, i.e. the orbital energy corresponding to the top of the 
valence band, is the exact ionization potential of the system [3]. This is known 
as the ionization potential theorem. Since adding an extra electron at the 
bottom of the conduction band would hardly change the density of a bulk 
system, the ionization potential theorem leads one to expect that the differ- 
ence in the Kohn-Sham orbital energies corresponding to the bottom of the 
conduction band and the top of the valence band should give the correct band 
gap of semiconductors and insulators. However it is well known that such an 
interpretation leads to gross underestimatation of the band gaps of semicon- 
ductors in comparison to the experiments. This is true even in solids like Ge, 
Si and the Ill-V semi-conductors which are not "strongly correlated". This 
underestimation of the band gap has been explained in terms of the derivative 
discontinuity [4,5] of the exchange-correlation energy functional. To circum- 
vent this difficulty, we propose to look at the calculation of the band gap as 
an excited-state problem. As such, the question we ask is: could we set up a 
DFT for excited states, calculate the total energy for appropriate excited-state 
and obtain the excitation energy as its difference from the ground state total 
energy. We may then go on to interpret the band gap in a semi-conductor as 
the difference of the energies between the ground state and an excited state 
where an electron in the highest occupied (HO) state of the ground state is 
removed from it and excited to the lowest unoccupied (LU) state. 

The mathematical basis of a density functional theory for total energies in 
excited states had been set up by Gorhng [6] and Levy and Nagy [7], based 
on the constrained search approach [8]. The theory has been put on a strong 
footing by Samal and Harbola [9]. The energy of the excited state of electrons 
in a solid may be expressed as a bi-functional of both the excited state density 
and the ground state density. 



where V{r) is the ion-electron potential. F[pex,Pgr] is a bi-functional of the 
excited state density Pex(i") and the ground-state density Pgr(i')- It is defined 
[7] from the constrained-search formulation : 



where the constrained search is carried out with the wavefunctions giving the 

excited-state density; T and V^e are the kinetic energy and electron-electron 
interaction potential operators respectively and the wavefunctions $ are or- 
thogonal to all the lower states determined by pgr. For example, in defin- 
ing F[pex, Pgr] for the first excited state, the $-s chosen would be such that 




(1) 



-^[Pex, Pgr] 



(2) 



2 



they are orthogonal to the ground state wave-function obtained from pgr. The 
exchange-correlation energy functional is then : 



-£^xc[Pex, Pgr] = -^[Pex, Pgr] " 7o[Pex] " 14i[Pex] 



(3) 



Here Tq represents the kinetic energy of a non-interacting system of electrons, 
suitably chosen so that its charge density is the same as that of the interacting 
system. The effective potential seen by the system of non-interacting electrons 
is variationally obtained and the single-particle wave-functions of the non- 
interacting system satisfies the Kohn-Sham equation : 



where V{r) — J2n '^('^ ~ f^n) are the positions of the ion-core labelled 

by n. The excited state density is obtained from : p{r) = J2i fi where 
the occupation number is or 1 depending on whether the state labelled 
by i is occupied or not. 

As in the ground-state theory, the exchange-correlation energy functional in 
excited-state theory too has to be approximated. The problem of finding an 
excited-state exchange-correlation energy functional is a challenging one and 
is further compounded by the fact that the functional is expected to be state 
dependent. Thus different functionals may have to be constructed for different 
classes of the excited-states. Recently, Samal and Harbola [10] have proposed 
a local-density approximation (LDA) for the exchange energy functional for 
a theory of excited-states with one gap. The functional has been applied to 
a large number of atomic excited-states to obtain very accurate excitation 
energies for single and double-electron excitations. The functional has precisely 
the form required to calculate band-gaps of semiconductors using excited-state 
DFT. The purpose of the present work is then to employ this excited-state 
DFT to obtain the band-gaps of a number of semiconducting and insulating 
materials. In the following we give a brief description of the methodology 
employed by us and then present our results. 

As we discussed earlier, the band gap in a semiconductor or an insulator can 
be thought of as the lowest energy required to excite one electron out of a 
many-electron ground state to one of the excited states of the system. Thus, 
in terms of the transition energy. 



where £^ex[Pex, Pgr] is the total excited-state energy functional of (1) minimized 
by varying pex and Eg^lpg^] is the ground-state energy corresponding to the 
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configurations mentioned above. 

To obtain tiie total energies we turn to the solution of the variationally ob- 
tained Kohn-Sham equation (4) [2]. The charge density is given by, 

oo oo 

p(r) = E/^ k^Wr with fi-^e 

1=1 i=l 

Here p is either the ground state or the excited state densities depending upon 
the occupation numbers /j. The total energy functional is : 

= Eks - En[p] + E^o (6) 

Here p is either or pcx depending on whether we are talking about the 
ground or the excited state; E^c is the "modified exchange" energy functional 
which includes the difference between the interacting and non-interacting ki- 
netic energies as well. This is of the form E^ciPgr] or £'xc[pex, Pgr] , again depend- 
ing upon whether we are talking about the ground or the excited state. E^i is 
the Hartrec energy and E^c is the contribution from the modified exchange- 
correlation energy functional given by the last two terms in the Eq. (6). Note 
that for the sake of notational clarity the bi-functional notation for the excited 
state functionals in the above equation has been dropped. 



2 The excited state exchange- correlation functional 

There have been many attempts at obtaining accurate but approximate forms 
for the exchange-correlation functional for the ground state. The simplest of 
them, the LDA, follows the Dirac [11] approach and proposes the functional 
on the basis of exchange energy for the homogeneous electron-gas. As men- 
tioned above, adopting the same approach, Samal and Harbola [10] proposed 
an excited-state functional for the exchange-energy and applied it to atoms 
to obtain accurate excited-state energies. In the following we discuss the con- 
struction of such a functional for semi-conductors and insulators and employ 
it to get the band gap of these materials as defined by Eq.(5) above. The 
accuracy or otherwise of our proposition will be evident from the applications 
to several semi-conductors and insulators. 

The standard procedure for setting up of the exchange-correlation functional 
is to start from the equivalent expression for the homogeneous electron gas 
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Fig. 1. The pictures of the ground state configuration of a HEG (left) and an excited 
state (right). 

(HEG). This is shown in Fig. 1. In the semi-conductor or insulator in its ground 
state, the electrons are filled up to the Fermi level, that is, in reciprocal space 
electrons occupy states labelled by wave vectors from /c = to kp. The same 
picture is true for the HEG, but in addition we have 

kl = 37rVgr(r) (7) 

where /5gr(r) is the almost uniform electron density in the ground state. 

On the other hand, for the lowest excited state in the semi-conductor or insu- 
lator, the electron in the HO state in the ground state configuration is removed 
from that state and now occupies the lowest unoccupied state. For the HEG 
this excited state configuration corresponds to the following : first, we shall 
call all states which were occupied in the ground state as well as in the ex- 
cited state : core states. If Pcore(i') is the electron density due to the electrons 
occupying the core states, then all states from k = to fci are occupied, 
where 

kl = 37r2peore(r) (8) 

Next, if Prem{'^) IS the electron density due to the electron in the HO state 
of ground state configuration which has been removed in the excited state 
configuration, then 

kl - kl = Stt^ p,,^(r) (9) 

Finally, if padd is the electron density due to the electron in the lowest excited 
state of the ground state configuration, which has been added in the excited 
state configuration, then 



- kl = Svr^ Padd(r) 



(10) 



All states between the labels ki and k2 which were occupied in the ground 
state become unoccupied in the excited state. Whereas all states between the 
labels /c2 and which were unoccupied in the ground state become occupied 
in the excited state. 



The ground state density is Pgr(i") 



Pc 



Pc 

(r) + Padd(r) 



+ 



Prcm ' 

r) - 



Pi 



and the excited state 
cm(r) + Padd(r). Note 



density is Pcx(r) 

that in an excited state of the system the electrons will occupy the /c-space 
differently from the ground state. 

As stated before, the standard starting point for building up the ground state 
exchange-correlation functional, whose exact form is in general unknown, is 
the exact expression for the HEG. We shall use the same idea for the excited 
state. This new functional will be the basis for our LDA generalized for excited 
states (MLDA). 

The exchange energy in the MLDA can be obtained by integration over the 
reciprocal space equi-energetic spherical shell surfaces. In a recent paper Samal 
and Harbola [10] have explicitly obtained the expression for the exchange 
energy functional for the excited state, as described above. We shall quote their 
results here and refer the interested reader to the above mentioned reference 
for the mathematical details. 

TpMLDA rpcore , padd i padd— core 

X X ' X ' X 



where 



^core 



^add 



E: 



add— core 



47^3 
1 



87r3 / 



dh 



"^(H - k^iks - k2) + {kl - kiyivL 
2{k3 - k2)kl + 2{kl - kl)ki + . . . 



{kl - klf In { ^^±^\ - (kl - klf J^'^^'^ 



k2 - kl 



M - ^1 



'h + k2 



It is easy to check from above that when ki = k2 = ^2 = kp, Eq. (11) reduces 
to the ground state exchange functional. Equations (8)-(ll) will then provide 
the basis for the calculation of the total energies. If we now go back to Eq. 
(6), we note that the expressions for the ground state and excited state total 
energies are : 



-C'gr--C/gr 
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£^ex = - Eu[PeA + E^'^'^iPe., Pgr] 



Here E refers to J2i fi^i where fi are the occupation numbers already defined 
and ej are the Kohn-Sham energies. Assuming no change in the orbital eigen- 
values corresponding to the ground- and excited-state configurations, which 
is a good assumption for bulk systems, the band gap will be given by the 
minimum value for different excitations of : 

AE^AE^' - AEu + {E^'''^[p..;p^] - E^^'^M} 

(12) 

where AE^^ is the conventional Kohn-Sham gap, A£^h is the difference in 

the Hartree energy corresponding to the excited-state and the ground-state 
densities, where both pcx and are calculated from the orbitals obtained 
from a single ground-state calculation. If the sum of the last two terms is 
small, we may expect the Kohn-Sham gap to be the band-gap. However, we 
expect the major correction to the conventional Kohn-Sham gap to come from 
these terms. What these terms accomplish is to replace the LDA exchange en- 
ergy corresponding to the excited-state density pcx by the MLDA exchange 
energy. The latter is tailor-made - and therefore more appropriate - for get- 
ting exchange energy for an excited-state of the kind considered for the band 
gap calculation. We note that the functional proposed by Samal and Harbola 
also has self-interaction correction (SIC) terms in it. For extended systems 
we expect these terms to be negligible and shall justify our assumption by 
application to a few of the systems studied by us. 



3 Application to Solids 

We shall base our calculations in solids on the tight-binding hnear muffin-tin 
orbitals method within the atomic-sphere approximation (TB-LMTO-ASA) 
[12]-[15]. In this method the muffin-tin spheres are inflated to the Wigner-Seitz 
cell volume and the interstitial contribution is neglected. The wave-function 
is expanded in a basis of linearized muffin-tin orbitals : 

l/jair, Ei^na) ^J2'^RLa{Ekna) xL(^' -^kna) YLif) 
RL 

Xf^(r, Ekna) = [ ^ma{r) + {E]^na " E^ri) ^Riair))] 
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Here, YL{f) are the spherical harmonics and ,^0- are the spinor wave-functions. 
The muffin-tin orbitals Xlo-(^5 E\^na) belong to an atom of type i?, the angular 
momentum L=(£, m), (pm<j{f) are the partial wave solutions in the muffin-tin 
spheres, ipma are the energy derivatives of the partial waves, both evaluated at 
the expansion energies Ei,iua used for linearization. Tail cancellation ensures 
that the muffin-tin orbitals labelled by R vanish outside the atomic spehers 
centered at R. The quantum states of the electrons in the solid are labelled 
by the quantum numbers K = (kncr). 

Prom the Kohn-Sham band-structure we can identify the HO and LU states in 
reciprocal space. Using these values of K we construct the removed and added 
electron densities. The core electron density is obtained from the atomic-like 
the core states. These expressions are then used in Equations (8)-(12) to obtain 
the band gap. We again point out that both the ground- and the excited- 
state densities are constructed from a single ground-state calculation. This is 
justified since one does not expect lowest unoccupied state to change much 
even if an excited-state KS calculation is performed. An advantage of this is 
that the computational demand of our calculations is pretty much at the same 
level as the ground-state calculations. 

Application of our methodology to molecules and clusters have already been 
carried out without any difficulty earlier [9]- [10]. However, in its application 
to sohds there have been two points of confusion. This is related to the fact 
that in a bulk sohd the extended Bloch-hke electronic wave-functions for a 
given quantum state are infinitely thinly spread over the whole solid. The 
charge densities corresponding to a single quantum state are also 0(1) while 
the charge density of the valence cloud is 0(N) [N is the number of atoms) ^ . 

There are two questions we must answer : First, are we really exciting a single 

electron out of a sea of N {N 00) electrons, from the HO to the LU state ? 
Secondly, can we use this infinitely small (0(1) as compared to 0(N)) change 
to produce the finite energy differences which we define to be the band gap in 
Eq. (12) ? 

In this section we shall argue that the way we obtain the energy differences 
within the LMTO technique by normalizing the solution within a Wigner-Seitz 
cell and considering only energies per cell addresses this problem successfully. 
We should also note that the arguments given is not confined to the LMTO, 
but holds good for any electronic structure method that are based on "cellular" 
techniques, e.g. KKR, APW or LAPW. 

As we have seen earlier, quantum states of electrons in a solid are labeled by 
K. If there are N {N ^ 00) atoms constituting the solid each contributing m 

^ A quantity is called 0{N") if it scales as A^" with the number N of atomic-spheres 
in the system. For example, total energy is 0{N) while the energy per atom is 0(1). 
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electrons to the valence electron cloud, then the total number mN of available 
states are occupied, each state by one electron. If p(r) is the total charge 
density of occupied electron states, then 

j d^r p(r) = mN integrated over the whole volume. 



In all electronic structure techniques which divide the solid into identical 
Wigner-Seitz cells (WSC) centered at the sites R such that : 

N 

0= y where f)i?f|^i?' = Vi? V ^ 

R=l 

where Q, is the volume of the solid and Q,r that of a WSC centered at R, the 
basis of representation are WSC centered functions of the form : 



xL(r) 



XLo-(^)^i('^)Co- for r e WSC centered at R 

(13) 

otherwise 



L is a composite index, which for the KKR or the APW and their linearized 
versions, LMTO and LAPW, denote the angular momenta indices (ima). 
These basis functions are centered at WSCs and do not overlap. The charge 
density also may be broken up as : 



N N 



P(r) = E E /k Pi(r) = E P^'i^) (14) 
?=i 

where p'^(r) = ' 



R=l K R=l 
^R 



p (r) for r e WSC centered at R 
otherwise 



Note that each quantum state labelled by K can be occupied by a single 
electron. Here /k is the occupation number of the quantum state labelled 
by K. For example, in the ground state mN of the lowest energy states are 
occupied, so mN of the /k are 1 the rest are 0. In a crystalline solid, the 
densities within each WSC labelled by R are identical. However, since these 
partial charge densities have the same structure as Eq. (13), and are outside 
the WSC labeled by R, we shaU retain this label to indicate this. Prom the 
above it is clear that, since p(r) integrates to mN over N WSC : 

J dhp^{r) = m and J dh p^{r) = 1/N (15) 



9 



In the output of the cellular methods (including the KKR, APW, LMTO,LAPV\' 
the charge densities per WSC are scaled such that they are normalized to 1 
within a cell. 



/ 



d'r p^(v) = 1 

So that the scaled charge densities per K labelled state which we get as output 
from the cellular methods are related to the true charge densities per state 
as : 

pS(r) = m)p^{r) (16) 



If we now remove or add a scaled charge density p^{r) integrating to 1 in 
a WSC, then that is equivalent to removing or adding a true charge density 
integrating to 1/N in that cell. The total charge density added or removed in 
the whole solid then integrates over the whole volume to 1. In other words, 
we are exciting 1 electron from the valence cloud in the sohd from the HO to 
LU state. 

The total energy expressions of Eq.(6) can also be broken up as : 

TV 



R=l 



We now go back to the way we construct the charge density of the excited 
state : We divide the total charge density of the ground state into a density 
of the 'core' which docs not change on excitation Pcorc(r) and that due to the 
HO state Prem(r)- We then remove Prem(r) and replace it by the density due 
to the LU Padd(r)- 

Pex(r) =Pcore(r) + Padd(r) = Pcore(r) + Prem(r) + (Padd(r) - Prem(r)) 

= Pgr(r) + (5p(r) 

We also note from Eq. (11) that when Sp ^ E^^^'^lpg,] = E^^^[pg,]. 
Combining these two, we get : 



AE[p] = £;^^^^[pex(r)]-i?"^^[pgr(r)] 

= E^^^^Mr) + 6p{r)] - E^^^[p,,{r)] 

N 

= E {<"^^[pfr(r) + Sp^ir)] - Ek^^ip^ir)]} 

R=l 
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+ ... 



,11^ 

J dh 6p\v) A[pl{v)\+o{{5p^{v)Y) 

Hr 

J dh 5p«(r) A[pJ(r)]+o({5p^(r)r) 



+ ... 



N 

E 



(17) 



Here o(a;") stands for terms which scale as a;™", m > n. Note that (5p'^(r) ~ 1/A^ 
so that all terms N o (^{6p^{r)}'^^ become negligible as iV — > oo. In this limit 
we may forget these latter terms. 

The last but one step follows from the fact that all WSCs are identical in a 
crystalline solid. We note from Eq.(15) that although the total energies on 
the right side of the above equation are individually 0{N), the difference AE 
should be ~ 0(1). Since 5p{r) scales hke 1/N and the summation of it! provides 
a factor of N, the functional 74[p^(r)] ~ 0{1). 

The energies calculated in the TB-LMTO are energies per unit cell. The energy 
difference calculated within the TB-LMTO is : 

AE\ft^] = / 5p^{v) A[pJ(r)] (18) 

^R 



The only way in which both A£'[p] and AE[p^] can be 0(1) is if the functional 
A[p\ has a specific scaling behaviour with p. Let A[p] = A[Np] — N°'A[p\, then 
q; = for the above to hold. In that case, it now follows that : 

AE[p^] = AE[p] = Eg 



The calculation in the solid is thus reduced to the calculation within a single 
WSC using the scaled charge density p(r) and the differences in energies per 
unit cell. We are exciting one electron from the entire valence cloud of the 
solid from the HO to LU state and the total energy difference is = 0(1). 
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4 Results and Discussion 



We shall first illustrate our calculations with a specific example, i.e. that of 
Si. The Fig. 2 shows the Kohn-Sham band structure of Si and the HO and LU 
states as obtained from the regular Kohn-Sham LDA calculation. In addition 
to these states, the calculations also give the Kohn-Sham gap AE^^. These 
states are then used to construct the removed and added electron densities for 
the excited state. When employed in Eq. (12), these densities yield the band 
gap of silicon. The numerical results are shown in Table 1 while a comparison 
is shown in Fig 3. 

15 
10 

I ' 
>. 
D) 

o 

lU 
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Fig. 2. The Kohn-Sham band structure of Si showing the HO and LU states used 
in our calculations 

As is clear from the Table 1 and Fig. 3, our method improves the Kohn-Sham 
band gap substantially, bringing it very close to the exact exchange [16] and 
experimental [22] band gaps. Interestingly the correction to the Kohn-Sham 
band gap is 0.52 eV and arises essentially from the exchange-energy difference 
between the ground- and the excited-state. This difference is very close to the 
derivative discontinuity of 0.58 eV in the exchange-correlation potential for Si 
obtained [17] from a GW calculation. 

Having obtained accurate results for Si, we have performed similar calculations 
for Ge, the III-V zinc-blende semiconductors GaP, GaAs, GaSb, InP, InAs and 
InSb, two zinc-blende oxides MgO and CaO, and one Wurzite structure oxide 
ZnO. The results are displayed in Table 1 and Fig. 3. It is evident from the 
results that our procedure yields uniformly accurate results for the band gaps 
of narrow band-gap to large band-gap materials. However, for semiconductors 
like InAs and InSb which actually turn out to be metals with zero Kohn-Sham 
gaps, our procedure is problematic. Therefore we did the initial calculation 
with a lattice constant slightly smaller than the equilibrium one, so that a very 
small Kohn-Sham gap appears and then redid our calculations. The estimated 
band gaps in the MLSDA are slight over-estimations. 

As mentioned before, since the band gap is calculated as an energy difference. 
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Table 1 

Comparison of the calculated band gap (AE) by the present work (LMTO-MLDA) 
with the TB-LMTO Kohn-Sham gap, the LMTO based Hartree-Fock (LMTO-HF) 
gap, the LMTO based exact exchange gap (LMTO-EXX) and the experimental 
(Exp) gap for diamond lattice Si and Ge, zinc-blende III-V semiconductors and 
wurzite oxides of Zn, Ca and Mg. 



System 
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LMTO 


LMTO 


LMTO 


Exp 






(LDA) 
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Si 
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0.49 


1.01 


5 6 [201 
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1.17 [22] 


Ge 
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0.53 
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0.74 [29] 
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1.62 


2.30 
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InP 
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0.71 


1.65 






1.42 [29] 


InAs 


6.04 [25] 


0.03 


0.61 






0.43 [29] 


InSb 


6.48 [25] 


0.01 


0.59 






0.23 [29] 


MgO 


4.21 [28] 


4.49 


7.64 


25.30 [21] 


7.77 [23] 


7.83 [24] 


CaO 


4.80 [28] 


3.35 


6.85 


15.80 [21] 


7.72 [23] 


7.09 [24] 


ZnO 


3.25,5.21 [27] 
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3.35 






3.44 [29] 
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Fig. 3. Band gaps calculated by different methods compared with the experimental 
band gaps of a series of semi-conductors and insulators. 
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we expect the difference in self-energy correction (SIC) between the added 
and the removed states to be small. We have carried out SIC corrections, as 
proposed by Samal and Harbola [9,10], for Si, GaSb and ZnO. For Si the gap 
changes from 1.01 eV to 1.02 eV; for GaSb from 0.94 eV to 1.06 eV and for 
ZnO from 3.35 eV to 3.37 eV. In all cases the correction is ~ 1% confirmig 
our assumption. 

Our correction is to the exchange energy functional alone. Corrections to the 
correlation part is being implemented in a subsequent communication. How- 
ever, even with the pure exchange correction major part of the correction to 
the Kohn-Sham gap is achieved. We are now applying this technique to a 
series of semi-conducting and insulating systems. 
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